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I. INTRODUCTION 

The problem of the coexistence of superconductivity 
and magnetism has attracted renewed interest during 
the last decade due to the discovery of unconventional 
and high-T c superconductors, in which such a coexistence 
seems to be realized. The presence of magnetic properties 
has been supported by experimental evidence in such ma- 
terials as high-T c ruthenocuprates, 1 Sr2Ru04, 2 ZrZn2, 3 
and UGe 2 ^ 

Experimental investigation of ferromagnetism coexist- 
ing with superconductivity is hampered by the fact that 
the spontaneous magnetic moment is screened by the 
Meissner currents in macroscopic samples. Therefore the 
usual methods, such as the Knight shift or muon spin 
relaxation, are only able to detect small remanent fields 
near inhomogeneities, defects etc. On the other hand, dy- 
namical measurements - e. g., excitation of spin waves 

- provide a direct probe of bulk magnetization and are 
sensitive not only to its magnitude, but also the direc- 
tion and degree of uniformity. There have been reports in 
the literature^ of microwave measurements in powder or 
ceramic polycrystalline samples made of materials with 
coexisting superconductivity and magnetism (SCFM). 

Theoretical research on macroscopic electromagnetic 
properties of SCFM's was initiated by Ng and Varma)Sil 
who have considered the bulk magnetic response and cal- 
culated the spin-wave modes for the spontaneous vortex 
phase. This regime, in which vortices are created in the 
sample without any applied magnetic field, is realized 
when the spontaneous magnetization is strong enough - 
namely, larger than Hd/Air, with H c \ being the first crit- 
ical field. However, in order to make direct comparison 
with microwave experiments, the finite-sample response 

- namely, the surface impedance - is needed. Unlike in 
usual metals, dielectrics, or paramagnetics, this quan- 
tity cannot be trivially related to the bulk properties in 
SCFM's, but requires a separate study. Hence the next 
step was done by the present author and Sonin^2ii£ who 
have considered the electromagnetic response of a thick 
slab, concentrating on the Meissner regime, in which 
static magnetic fields are screened in the bulk. In this 
regime, static magnetic measurements are difficult, and 



hence microwave experiments may be of particular inter- 
est. Spin-wave excitation and propagation were consid- 
ered for a SCFM slab with the magnetization perpen- 
dicular to the surface (perpendicular geometry). 8 While 
the analysis was initially carried for the case of Landau- 
Lifshitz (LL) SCFM's with spin magnetism, the results 
have been generalized to the case of p-wave superconduc- 
tors with orbital magnetism.? It was shown that these 
materials are in general described by more complicated 
magnetic dynamics than LL dynamics. However, two 
specific cases - namely, when the waves propagate either 
parallel or perpendicular to the spontaneous magneti- 
zation direction - can, in fact, be described by LL dy- 
namics, though corresponding to different values of the 
magnetization, thus demonstrating the relevance of the 
analysis to unconventional superconductors with orbital 
magnetism. Also, it was found that SCFM's can sup- 
port a surface wave^S and its spectrum was analyzed as 
compared to normal ferromagnets. 

In the present work, I extend the microwave response 
analysis started in Ref. |g, including both the perpen- 
dicular and parallel geometries (when the magnetization 
is perpendicular or parallel to the surface). Noteworthy, 
the response is shown to exhibit a resonance related to 
excitation of the surface wave. The influence of applied 
magnetic fields is also considered, leading to completely 
different effects in the perpendicular and parallel cases. 
In the former case, the field penetrates the sample, push- 
ing it into the vortex state. In the latter case, on the other 
hand, the field is screened inside the Meissner layer and 
the system loses its uniformity. 



II. THE MODEL 

We consider a SCFM in form of a slab which is thick 
enough so that it is possible to neglect reflection of waves 
from the second boundary. The material is assumed to 
possess a ferromagnetic anisotropy of the easy axis type, 
with the direction of the easy axis being chosen as the z 
direction, so that the equilibrium magnetization is Mo = 
Mq z. For AttMq < H c i, which we assume, the material is 
in the Meissner phase, except when an external magnetic 
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field is applied perpendicular to the sample surface (in 
which case the material is in the vortex state). The free 
energy functional for the material is given by£*& 
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where M is the magnetization, M. x ,y is the magnetiza- 
tion component normal to the easy axis, di — d/dxi, 
a > is the anisotropy parameter (note that it has a dif- 
ferent normalization here than in Ref. 8), (f) is the phase 
of the superconducting order parameter, <3?o = hc/2e is 
the magnetic flux quantum, A is the vector potential, 
and B = V x A is the magnetic induction. The ap- 
plied static magnetic field is parallel to z and so is the 
static component of the magnetic induction Bo- The pa- 
rameter 7 characterizes the stiffness of the spin system. 
The anisotropy is assumed to be strong enough - namely, 
da = a — 4-7T > - so that the system is stable (or at 
least metastable) against flips of the equilibrium mag- 
netization Mo. In accordance with the micromagnetism 
approach, the absolute value of M is assumed to be con- 
stant, and hence terms dependent on |M| were omitted. 
We consider a single-domain configuration since this con- 
figuration is either stableii (for strong enough supercon- 
ductivity 2tt\ < 7a 1 / 2 ) or at least metastable (in the 
opposite case) with respect to domain formation, due to 
the anisotropy, as explained above. In our treatment, 
the dissipation initially is not taken into account explic- 
itly, so its only manifestation is the absence of reflected 
waves from the other surface of the sample. Then, having 
obtained the result, we generalize it, including explicitly 
Ohmic normal currents. 

The spin dynamics is governed by the Landau-Lifshitz 
equatioi 
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where g is the gyromagnetic factor. Since we are con- 
cerned here with motion near the equilibrium, we can de- 
compose M = Mo 4- m and B = Bo + b, where m _L Mo 
and b are dynamical parts. For a given frequency Eq. P)l 
yields 

-ium= -</M x{(a-7 2 V 2 )m-b}-gB xm. (3) 

Another equation, needed for the dynamical part of 
the magnetic induction b, is the London equation 



(1 - V 2 A 2 ) b = 4ttA 2 V x V x m 



(4) 



The same equation also determines the decay of the static 
component of the magnetic induction parallel to the sur- 
face inside the sample: 
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FIG. 1: The setup in the perpendicular geometry. 



In the regions where Bo is significant, translational sym- 
metry is broken, so the resulting equations are quite com- 
plicated. However, if for some reason Bo = (e.g. in the 
bulk of the sample, where Bo has already decayed, or by 
application of an appropriate external field) , the transla- 
tional invariance is restored, and the solutions arc plane 
waves ~ e «(k-x-o)t)_ Then equation of motion for the 
magnetization is 
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where VL = to/ gMo is the frequency in magnetic units and 
mi is the divergence-free component of m. Note that re- 
tardation effects can be safely neglected here, since the 
displacement currents are negligible in comparison with 
the superconducting currents as long as the frequency is 
much smaller than the plasma frequency, which we, of 
course, assume. This is in contrast to insulating fer- 
romagnets, where mixing of spin and electromagnetic 
modes can be important. To proceed further, we will 
consider two geometries, in which the easy magnetization 
axis is perpendicular and parallel to the sample surface. 



III. PERPENDICULAR GEOMETRY 

Here we consider a case when the easy axis is perpen- 
dicular to the sample surface, shown in Fig. ^ I n this 
geometry all quantities vary along the z axis, so m is 
divergenceless: mi = m. At zero applied magnetic field 
the static induction Bo vanishes and the system is trans- 
lationally invariant. We first consider this simpler case 
and then generalize the results for the case of nonzero 
applied magnetic field. 



A. Zero applied field 

In this regime, Bo = and m± = m, so from Eq. © 
the spectrum of spin waves is 



n = ± a + 7 2 fc 2 - 
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FIG. 2: The spin-wave spectrum in the perpendicular geome- 
try in two regimes: stiff, e = -\/47r + 0.2, a = 4tt + 50 (dashed 
line) and soft, e = V^tt - 2.35, a = 4tt + 0.013 (solid line). 



The spin-wave modes are circularly polarized, and the 
two signs above correspond to two senses of the polar- 
ization: = rrv- )(x =F iy). The polarization is de- 
termined by the incident electromagnetic radiation due 
to continuity of the fields across the slab surface. Only 
positively polarized waves can propagate inside the sam- 
ple, and further on we mainly focus on this polarization. 
Depending on the ratio e = 7/A, the spectrum can take 
two different forms, as shown in Fig. if e > V^tt (high 
stiffness, strong superconductivity), the spectrum has a 
minimum (threshold for spin wave propagation) Qf = a 
situated at k = 0. In the opposite case, the minimum fre- 
quency is at a finite wave vector k m = X~ 1 (^/4Tr/e — l) 1 / 2 
and has a lower value O m = a — (\/47r — e) 2 . In this 
regime, waves with wave vectors satisfying \k\ < k m have 
negative group velocity duj{k)/dk. In order to carry the 
energy away from the boundary, these waves should have 
wave vectors directed toward the boundary. 

For given frequency and polarization, two spin-wave 
modes with different wave vectors are excited by the in- 
cident radiation. For each mode, the corresponding elec- 
tric and magnetic fields inside the sample are found from 
the Maxwell equation dh/dt = — cV x e and the London 
equation (0J: 
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The total fields are given by superposition of different 
modes. In order to find a proper superposition of the two 
spin-wave modes inside the sample, an extra boundary 
condition is needed, in addition to the usual continuity 
of tangential components of the electromagnetic field at 
the sample surface. This additional condition should be 
imposed on the magnetization. The simplest possible 
boundary condition 12 is dm(z)/dz — on the surface, 



FIG. 3: Re<^ + ' for the perpendicular geometry in the stiff 
regime, e = + 0.7 (dashed line) and soft regime, e = 



4ir — 1.7 (solid line). In both cases, a — 4tt + 0.2. 



which physically means absence of spin currents through 
the sample surface. In the Fourier space this condition 
gives 



feimi + k 2 m 2 = 



(9) 



This equation together with Eq. (JHJ formulates the 
boundary problem, from which all fields in the slab can 
be found. 

Now we turn to calculation of the microwave re- 
sponse, which is quantitatively expressed by the surface 
impedance and whose real part is proportional to the en- 
ergy absorption by the system^ For circularly polarized 
fields, it is given by 
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Using Eqs. 
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and 10, we obtain, after simplification 
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where qi = kiX. Then, using the dispersion relation, 
Eq. (J7J, one can solve for the wave vectors and obtain a 
general expression for the surface impedance as a func- 
tion of the frequency: 
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where Q = i\/±Q — a. It can be observed that 
pure imaginary since spin waves with negative circular 
polarization cannot propagate, and it shows regular be- 
havior without special features. On the other hand, for 
the positive polarization, the impedance has a real part 
satisfying the inequality Re(^ + ) > 0, as it should, and 
similar to the spectrum, it shows a qualitatively differ- 
ent behavior in the two regimes depending on the value 
of e as can be seen in Fig. Most importantly, it has 
several singular points (thresholds) related to spin-wave 
propagation properties, which can be especially useful for 
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comparison with experiments. Below we list these points 
and analyze the behavior of (^ + ) around them. 

(i) The case e > V^n (stiff spin system). Then fi(fc) 
is a growing function and the threshold frequency for 
spin- wave propagation is Qf — a. This frequency cor- 
responds to uniform oscillations in SCFM, at which the 
electric field vanishes inside the sample. Near Clf one 
wave vector is small {k\ cx ^JQ, — Oj), and then the sur- 
face impedance can be simplified: 
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Thus, ReC (+) is zero below the threshold f2/ and has a 
square-root behavior above it 

(ii) The case e < \J\~k (soft spin system). In this 
regime the threshold frequency has a lower value fi m = 
a — {\[A~k — e) 2 < flf. At this frequency the electric 
field on the sample surface vanishes (but unlike at ttf , it 
is nonzero inside the sample). Near O m , the two wave 
vectors are close in absolute value to k m , but oppositely 
directed: k\ + ki oc y/Q — f2 m . Then the impedance is 
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Thus, ReC'" 1 "- 1 has a square-root singularity at Q m . 

In addition, there is still a singularity at fi/. The be- 
havior there, however, is different than in the stiff case: 
namely, Re(^ + ) [ s nonzero both above and below fi/. In 
order to obtain it, Eq. (|13l) is not sufficient, and one has 
to go to the next order of smallness. Then 
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so for a soft spin system, Re£( + ) vanishes at f2/, showing 
a square-root behavior below it and a linear one above 
it. 

In addition to these points, where the impedance van- 
ishes, there is a frequency fifc = da — e 2 /2 + ey / An + e 2 /4 
at which <^ + ) has a pole for both stiff and soft regimes: 
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Then Re has a 5-function peak at ujh which acquires 
a finite width if small dissipation is added to the system. 
This peak corresponds to the ferromagnetic resonance 
for a SCFM slab, and it is also related to the surface 
wave supported by the system, whose branch starts at the 
same frequency lu^ in the limit of infinite light velocity. 10 
Indeed, the oscillation excited at this frequency by an 
incident radiation is identical to the oscillation produced 
in the slab by the long- wave surface wave. Finally, at 
large frequencies f2 >> a the magnetic properties stop 
playing a role and the surface impedance becomes just 



the impedance of a plain superconductor, £ = —iuj\/c, 
with the real part being a small correction ~ cj~ 3 / 2 . 

Our analysis can be easily generalized to take into ac- 
count dissipation due to normal currents. For this, nor- 
mal current J n = <r n E should be added to the supercon- 
ducting current in the Maxwell equations. Here a n is the 
normal conductivity and the frequencies are assumed to 
be not very high, so the dispersion of the conductivity 
can be neglected. As a result, in the London equation 
(|4*|). a renormalized complex X n should be used instead 
of the London penetration depths 
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where <5 = c/ \J2ii(j n u) is the skin depth. Consequently, 
the spin-wave spectrum and the surface impedance are 
obtained from Eqs. (J7|, (|11|) . and l|12|) by analytic con- 
tinuation using complex A„ and e. Then, in general, the 
distinction between the soft and stiff cases disappears, as 
Re £( + ) loses its singularities and becomes nonzero at all 
frequencies, except at ujf, where it shows a square- root 
behavior both above and below ujf. However, if the nor- 
mal currents are small, as happens in a superconductor, ( 
behaves approximately in the same way as in the nondis- 
sipative situation discussed above, with the singularities 
smeared by a small amount ~ A 2 i5 -2 . In particular, the 
peak at ujh acquires a finite width and can be observed. 

The opposite limit A — * oo and A 2 = z<5 2 /2 corre- 
sponds to a metallic ferromagnet. In this case Re£ again 
shows no features (except vanishing at Wf). If, how- 
ever, the conductivity is not very high, 7 << 6, then 
the frequencies lu^ and uj m (calculated with a complex 
e = 7/A„) merge and the impedance shows a sharp 
peak there (with width ~ j/S). This property has been 
widely used in ferromagnetic resonance measurements 
111 metals»i^ Thus, superconductivity modifies the mi- 
crowave response in a qualitative way, which can be used 
for detection of the superconducting transition. 

Finally, we would like to make the following comment. 
By analogy with metals, one might think that the limit 
A — * 00 corresponds to an insulating ferromagnet. This 
is not so, since in the description of the latter retardation 
effects are important and should be taken into account. 
Accordingly, the correct limit describing a magnetic in- 
sulator is A~ 2 — ► —(uj/c) 2 . 



B. Finite applied field 

When a finite magnetic field Bq << H C 2 is applied 
perpendicular to the sample surface, it penetrates the 
sample, where the vortex phase is formed. In the pres- 
ence of vortices the London equation is modified since 
the equation for the superconducting current density j s 
now becomes 
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where u is the local displacement of the vortices from the 
equilibrium position. To complete the description of the 
problem in this regime, one needs also an equation for 
the vortex dynamics. The simplest reasonable model is 
that the vortex motion is governed by a balance of the 
Lorentz force and a viscous drag. In the presence of nor- 
mal currents, one should be careful about the form of 
the Lorentz force, so that the Onsager relations are not 
violated. 16 A plausible assumption is that the Lorentz 
force depends only on the superconducting current den- 
sity. With all these considerations taken into account, 
the vortex dynamics is described by 
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where r\ is the viscous drag coefficient. Substituting this 
equation in the previous one, we obtain 
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where A 2 = A 2 + iBq^q/Attuit]. This equation formally 
coincides with the London equation in the Meissner state, 
so in our approximations the presence of vortices amounts 
to a replacement of the London penetration depth A by 
the complex A„. Then the presence of normal currents 
is taken into account by a further replacement 16 A~ 2 — > 
^nv — \ 2 ~ 2«<5 2 . 

In addition, it is evident from Eq. @ that instead 
of the bare anisotropy parameter a a modified a — 
a + Bq/AttMq should be used. Thus, the effect of the ap- 
plied field is: (i) to shift the position of singularities in the 
response (this can be useful since it allows to take mea- 
surements at constant frequency by scanning the applied 
field) and (ii) to add dissipation to the system, smearing 
the singularities and eventually pushing the response into 
the metallic regime. 



IV. PARALLEL GEOMETRY 

Now we consider a situation when the easy axis is par- 
allel to the sample surface, as shown in Fig.0J The direc- 
tion perpendicular to the surface is chosen to be x, along 
which all quantities vary, and then = ym y . Again 
we start from a simpler case when the magnetic induc- 
tion vanishes inside the sample, Bq = 0. After that, we 
consider a more complicated case when there is a nonzero 
induction near the sample surface, decaying in the bulk. 



A. Zero static induction 

The regime in which the static magnetic induction van- 
ishes inside the sample can be realized by application 
of an external field Ho = — 47rMo. Strictly speaking, 
this configuration is not stable, since flipping Mp in the 
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FIG. 4: The setup in the parallel geometry. 



opposite direction would lower the energy. However, it 
is metastable provided the anisotropy is strong enough, 
a > 47r. From Eq. JBJ the spectrum of spin waves now is 
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It can be seen from this relation that in the parallel ge- 
ometry an incident electromagnetic wave excites in the 
sample three spin-wave modes. Since the dynamical part 
of the magnetization has only one component in the plane 
of the sample surface, m y , the spin wave in this geometry 
is excited by a linearly polarized incident radiation with 
the magnetic component in the y direction. The spec- 
trum has the same qualitative features as in the perpen- 
dicular geometry, although the specific values of different 
parameters are, of course, different. Thus, the transition 
between the soft and stiff regimes occurs now at e = y/2n, 
while the wave vector k m || and the frequency uo m || can be 
found by minimizing Eq. (|21|) , which requires solution of 
a cubic equation. Note, however, that the frequency of 
uniform oscillations Q f remains equal a in this geometry. 

The boundary condition, Eq. 10, now gives two inde- 
pendent equations for x and y components of the magne- 
tization. These two components are related by Eq. ©: 



iSlrriy = (a + r y 2 k 2 )m a 



(22) 



Using this, the boundary conditions can be written as 



kj m x j = 
^ kf m xi = 0, 



(23) 



where the summation is over the three spin-wave modes. 
These conditions determine the relative amplitudes of the 
modes and together with the dispersion relation specify 
the problem completely. 

The surface impedance for this case of linearly polar- 
ized incident wave is given by 
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FIG. 5: Re£ for the parallel geometry in the stiff regime, 
e = \/2ty + 0.1 (dashed line) and soft regime, e = — 1. 
(solid line). In both cases, a = 4n + 0.3. 



Then the calculation goes along the same steps as for the 
perpendicular geometry, though it is more complicated 
due to the presence of three modes. Using Eqs. ©, JUJ, 
and J23Jl, one obtains, after some algebra, 
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where cij are symmetric combinations of the normalized 
wave vectors for different modes (and, as before, % = 
ki\): 



ai = 9l +92 + 93, 
a-2 =919293, 

= 9i 92 + 9293 + 939i- 



(26) 



These combinations can be found from the dispersion re- 
lation, Eq. H21f> . thus providing the dependence of ( on 
the frequency. As can be seen from Fig. EI Re C shows the 
same qualitative behavior as in the perpendicular geom- 
etry around its singularities. The same is also true about 
the high-frequency behavior of £, given by — iloX/c with a 
small real correction ~ w~ 3 / 2 . Below we give expressions 
for C, near the propagation thresholds. 

(i) In the stiff regime e > \/2tt, the impedance near the 
threshold Qf = a shows a square- root singularity 
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(ii) In the soft regime e < V2tt, the threshold is at 
u m ||, again with a square-root behavior around it: 
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Finally, Eq. I|25|l has a pole at a frequency || , at which 
ReC shows a peak, just like in the perpendicular geome- 
try. 



B. Finite static induction 

If the applied magnetic field Ho (parallel to the surface) 
is not exactly equal to — 47rMo, then there is a nonzero 
magnetic induction penetrating the sample, 



B (a;) = (4^M + H )e- 
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In this case the translational invariance of the system 
is broken, so the spin-wave modes are not plane waves 
anymore. This considerably complicates the problem, 
since these modes are given, instead of Eq. J3TJ, by a 
sixth-order differential equation 
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where / = (Att + Hq/Mo) is the magnetic induction at the 
sample surface divided by Mo. Then the general solution 
of the problem including finding the modes, matching 
them by the boundary conditions, and calculating the 
surface impedance would be very difficult. Instead of this 
we limited our attention to two extreme regimes in which 
the problem can be simplified. Namely, we considered the 
regime of very low stiffness (very soft regime), e << 1, 
and the opposite very stiff regime e >> 1. In the former 
case, the Meissner layer is wide and the static magnetic 
field there is very smooth, and therefore the quasiclassical 
approximation can be applied; in the latter, the Meissner 
layer is extremely thin and its influence is weak. 



1. Very soft regime 

In this regime e « 1, and the dispersion of free spin- 
wave modes (those outside of the Meissner layer) can be 
obtained by simplification of Eq. (|21|) : 



fl 2 = (a + j 2 k 2 ){Sa + j 2 k 2 ) + 
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Then the threshold frequency f2 m || for spin-wave pene- 
tration is given by 



n 2 m u = aSa + 4 ev /27ra(a - 2tt) + 0(e 2 ) (32) 
and the corresponding wave vector is 
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Of the three spin-wave modes corresponding to a general 
frequency f2, two have short wavelength ~ 7, while one 
has long wavelength A. However, for frequencies near 
Q m || , the picture changes: namely, one has a short-wave 
decaying mode with k 2 = — 2 (a — 2ir) r y~ 2 and two modes 
with intermediate wavelengths ~ VtA- Since in this case 
all wavelengths are much shorter than the scale of varia- 
tion of the "potential" [provided by the field Bo (2;)], it is 
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possible to treat the spin- wave equation i|30|) by the qua- 
siclassical approximation (its validity is discussed below 
in more detail). 

Thus we will be interested in finding £ (and in par- 
ticular, its real part) for frequencies slightly above the 
threshold, < Cl — Cl m \\ ~ e and in the presence of weak 
magnetic induction Bq ~ e-Mp directed either parallel or 
antiparallel to Mo. We start by discussing ( in the ab- 
sence of Bq. From the free- wave spectrum, Eq. I|31|) . one 
finds relations between the wave vectors k\ and k 2 : 



kl 



k\ 



k 1 k 2 — 



Cl 2 — a5a 
2(a- 2tt) 
2ita 



7 



2tt 



7 



4 A- 



(34) 



while the short-wave-mode wave vector k% has been spec- 
ified earlier. From these relations it is not difficult to find 
the wave vectors and substitute them into Eq. 1|25[) . thus 
obtaining 



c 



ci 2 



n l ii 



(35) 



where the resonance frequency Cl^ || is given (in the very 
soft regime) by 



Cl\ n = aSa + 2e^2ixa(a - 2tt). 
We recall that Eq. I|35|l is valid only near 



Eqs. hold, 
reduces to 

c 



To lowest order in ui — 



J m ||' 



(36) 
where 



|| 



Eq. 



U m || A y/ Cl m || (Cl - Cl m || ) 

c €\fa. — 2n 



(37) 



which coincides with Eq. (|28|l in the very soft (e << 1) 
limit. 

Let us now turn on the magnetic induction Bo (a;); 
then, its presence can be effectively described by an in- 
dependent anisotropy factor a(x) = a + Bq(x)/Mq. Ac- 
cording to the quasiclassical approximation, the corre- 
sponding spin- wave modes are found from the free modes 
by using x-dependent wave vectors. Note that for the 
short-wave mode 3 this dependence can be neglected 
since it is small and regular. On the other hand, for 
the modes 1 and 2, this dependence is crucial, since it 
determines properties of the wave propagation such as 
tunnelling. Indeed, substituting a(x) into Eq. (|32|l . one 
obtains the cc-dependent threshold frequency 



^||(z) = ^||+2(a-2^)/ e ^/\ 



(38) 



which can be considered as an effective potential as shown 
in Fig.|Sl In the region where CI > Cl m \\(x), modes 1 and 
2 are propagating, while if CI < Cl m »(x), these modes 
can only tunnel under the potential barrier. The point 
Xt between these two regions is a classical turning, or 
reflection point. It is defined by the condition 



CI = Cl m || (x T 



(39) 




Tunneling Turning Propagating 
. point modes 



FIG. 6: Coordinate-dependent threshold frequency u> m « (x) 
serving as an effective potential for Bo > (solid line) and 
Bo < (dashed line). 



The field-free threshold f2 m || = f2 m ||(x = oo) remains 
the propagation threshold so that only radiation with 
Cl > Cl m || can penetrate (here we ignore the peak at 
u>h\\)- When B(x) < 0, this penetration is classically 
possible, since Cl > Cl m \\(x), while when B(x) is positive 
and not too small, so that Cl m ||(0) > Cl, it is realized by 
tunneling. 

The validity of the quasiclassical approximation is of 
special concern in this problem. Indeed, as was men- 
tioned above, one can hope to use the quasiclassics only 
near the threshold Cl m » . On the other hand, Cl m u (x) = Cl 
means that the motion is near the turning point, around 
which the quasiclassics breaks down. The usual condition 
of validity d x k(x) << k 2 , is not sufficient here, since, 
unlike in quantum mechanics, the turning point corre- 
sponds not to zero but a finite wave vector k m so that 
Cl(k) — Cl m || oc (k — k m ) 2 (for the mode with Re A; > 0). 
Instead, a more stringent condition is required, 



d5k 
dx 



« 5k 2 



(40) 



where 8k = k — k m . There is also another dangerous 
frequency given by Cl 2 _ — aSa — 4e-y/27ra(a — 2n), at 
which k 2 = — fc^, the group velocity vanishes and the 
quasiclassics breaks down again (this frequency has no 
quantum-mechanical counterpart). However, it is situ- 
ated deep inside the classically forbidden region, and we 
ignore it, assuming Cl to be always larger than it. Then 
the above condition reduces to 



efe 



c/X 



« 



n 2 



Cl 2 Jx) 



a 



2tt 



3/2 



(41) 



(i) For B(x) < 0, the effective potential bends down. 
Then, as long as the quasiclassics holds, the wave prop- 
agation inside is described by two outgoing modes 1 
and 2 (plus the evanescent short-wave mode 3) with IE- 
dependent wave vectors found from Eq. I|34|l . Then £ can 
be calculated exactly as in the field-free case, but with a 
local value of a at the surface, a(x — 0) = a + f. The 
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result is 



'•>m ||- 



n 2 



n 2 „ 

m || 



2(a-27r)/ 



ft 2 



u h ii 



2(«-2tt)/ 



(42) 



which has just the field-free form shifted to the left 
(/ < 0). One might worry about this result that it gives 
real Q for f2 < tt m || . In fact, as £1 — > f2 m || , the quasiclas- 
sics breaks down at some depth inside the slab, and this 
result loses its validity. This is expressed in the appear- 
ance of significant over-barrier reflection. Then the cor- 
rect frequency dependence becomes Re£ oc y^fi — Q m || . 
However, as long as il is not too close to J7 m || [e.g., for 
/ ~ fl — fi m || ~ e, when condition (|41l) is satisfied for 
any x], the result given above is valid. 

(i) For B(x) > 0, the effective potential bends up, 
forming a barrier. For frequencies f2 m i| < 51 < fi m |i(0), 
the wave penetration is realized by tunneling of the 
modes 1 and 2, which means that each of these modes 
consists of a dominant decreasing part and a small (sub- 
dominant) increasing part (while the mode 3 remains un- 
affected). The tunneling occurs until the turning point 
xt, after which the modes become propagating. In order 
to calculate Im £, one may neglect the tunneling, thus 
assuming that all the energy is reflected from the bar- 
rier. Then only the dominant part of the modes should 
be retained, the problem becomes again formally anal- 
ogous to the free case, and Im£ is given by Eq. I|42|) 
with / > 0. Calculation of Re£ is more difficult, since 
it requires higher precision. It can be done directly by 
keeping both the dominant and subdominant parts and 
solving a boundary- value problem with five components. 
However, it is easier to solve the problem using the fact 
that ReC is proportional to the energy flux: 



ReC = — P\hy(P)\~ 



(43) 



The flux may be calculated in the region after the turning 
point, where it is given by 



i=l,2 



Um\m xi (x)\ 2 dil 2 (k[x] 



4a 



Ok 



(44) 



with Q,(k) given by Eq. (|31Jl . The transmitted ampli- 
tudes may be related to the amplitudes at the boundary 
by a procedure analogous to the quasiclassical calcula- 
tion of transmission coefficient in quantum mechanics, 17 
neglecting the subdominant parts. When this is done, 
one has to relate different mode amplitudes at x — us- 
ing the boundary conditions, Eq. (|23|) . Substituting the 
results into Eqs. and (3^ , one obtains 



ReC 



8Aw. 



m || 



e 3 / 2 (a - 2^) 3 / 4 (2™) 5 / 4 y / ft 2 i|| (0)-r! 2 



[n 2 - n 2 h u (o)] 2 ^n 2 -ni(o) 



(45) 




lOmu Whii (0) 



Wa,, (0) Cd 



FIG. 7: Effect of magnetic field in the parallel geometry, very 
soft regime. Re£ for / = (thick line), / = 0.05 (thin solid 
line, its scale given to the right), and / = —0.05 (dashed line), 
which is extrapolated to have a square-root behavior near ui m 
(dotted line). Note that the values of ui h ii (0) and io m n (0) are 
given for a positive field (/ = 0.05). Parameter values for the 
plots are a — 4tt + 1 and e = 0.01. 



where for all threshold frequencies, SI 2 (0) = tt 2 + 2(a 
2tt)/ and 



t = exp 



-2 / Imfc(ir) dx 
Jo 



(46) 



is the tunneling exponent. Using Eqs. I|34|l and l|39|) . one 

finds 



t = exp 
with 



-2e"V/e" iT (y^r - 1 - arctan y/e& - 1 

(47) 



2/ 



2tt 



n 2 - n 2 „ 

m || 



(48) 



Note that £ may have a peak of finite width if the field 
is strong enough so that io h || (0) > oj m || . The peak is cut 
off at Q 2 — S! 2 || (0) ~ e 3 / 2 i, when the subdominant part 
becomes significant and the tunneling is no longer small, 
as was assumed. On the other hand, the singularity at 
fi_(0) is disregarded, since, as we mentioned above, this 
frequency is assumed to be always smaller than £l m || . 
The calculation remains valid for il going down to H m || . 
The reason for this is that unlike in the previous case, the 
quasiclassics now does not need to hold everywhere in the 
sample (it is anyway broken around the turning point), 
but only in some finite region under the barrier and deep 
inside after xt- It can be seen from Eq. (|41|l that these 
requirements are satisfied for < Q — fl m m < / ~ e. 
Plots of ReC for Bq = 0, Bq < 0, and Bo > are shown 
on Fig. [7| One can see that for the antiparallel field 
(B < 0) the response is not affected significantly, unlike 
for the parallel field. 
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2. Very stiff regime 



and h,. 



In this regime e >> 1, which corresponds to the limit of 
very strong superconductivity. Thus the Meissner layer 
is very thin, and its influence amounts to modification of 
the boundary conditions. The threshold for wave propa- 
gation in this regime is Qt — a. Again we start by con- 
sidering the field- free situation first. The free spin- wave 
spectrum, Eq. H21|) . yields in this regime one short-wave 
and two long- wave modes: 



-A" 



k 2 

K 2,3 



(49) 



Substituting them into Eqs. (|25|) and J2HJ), one finds the 
field-free result 



c 



c 



2tt 



(50) 



Thus to the leading order £ is pure imaginary and corre- 
sponds to the surface impedance of a usual nonmagnetic 
superconductor, as it should in this regime of strong su- 
perconductivity. Note that this expression is not valid 
very near the threshold - namely, for Q, — a < (T 2 - and 
hence it does not capture the special behavior there, such 
as the square-root singularity £ ~ — a and a pole at 
Cl h || (in this regime uj h || and ujf are close to each other). 

Now we introduce the magnetic field, Eq. (|29[1 . which is 
not assumed small, and calculate the change in Re £ due 
to it. As will be shown below, the change comes only as 
a small correction ~ e~ 2 to the field-free results; hence, 
small terms ~ e~ 2 should be retained in the calculation. 
To start with, it is useful to rewrite the spin-wave equa- 
tion, Eq. i|30|) . in terms of a rescaled coordinate £ = x/X: 



n 2 (l-dl)m x = (l-df)(a + /e 



2q2n2 



rn. 



+4ird|(a + /e _t - e z df)m x . (51) 

This equation contains a large parameter e, which makes 
it possible to find the solutions as an expansion in e . 
As in the field-free case, one has a short-wave mode 



(52) 



and two long- wave modes 

m x 2,3(0 ~ m 2,3 exp(«Q2,30( 1 + / e ~ 2 e~ C ), (53) 

where we have only retained corrections related to the 
field (~ /). Similarly one finds the modes for m y , 

< / fe- 2 e-Z 



m yl (0 



2 m\e 



1 



m y 2,3(0 = =Rm>,3 exp(iq 2 , 3 £,)( l + ^ 2 / e ^). (54) 

where in the second equation the upper sign is used for 
mode 2 and lower for mode 3, and analogously for b y , 



4 mie 



in ' 

by 2,3(0 = T47rim 2 , 3 exp(iq 2 , 3 £) [ q\ 3 



, e- 2 fe-Z 

2*92,3 



(55) 



Km 



4 -£ 



h y 2,3(0 = ±47riTO 2 ,3exp(ig 2 ,30 1 



2iq. 



2.3 



.(56) 



Now one should use the boundary conditions in order 
to relate the amplitudes of different modes. It would 
be incorrect to use the conditions Vm(0) directly since 
differentiation would reduce the relative precision in the 
long-wave modes and the result would only have a preci- 
sion of e . Instead, one should integrate the equations 
for magnetization components, Eq. over the Meiss- 
ner layer and use the above condition in the integral of 
V 2 m, leading to 



[fe ? m y (£) - b y (£)] dx = e 2 (q 2 m 2 - q 3 m 3 ), 



fe e m x (£) dx = e z (iq 2 m2 + iq3m 3 ). 



(57) 



Substituting the expressions for the magnetization com- 
ponents and b y , one obtains relations between the mode 
amplitudes: 



n 



mi 1 



e- 2 f 



3=2,3 



±™ 3 ( -iQj + e~ 2 /(l + iQj) + 2tt 



iq 3 



m 2 (iq2 - e 2 f) + m 3 (iq 3 - e 2 f) = 0, 



(58) 



where again the upper and lower signs in the first equa- 
tion refer to modes 2 and 3, respectively. These relations 
are used in an expression for the energy flux in the region 
deeper than the Meissner layer, 



P = Lu\e 2 \m 2 \ 2 Req2 



(59) 



(where only mode 2 contributes as it is the only propa- 
gating mode), and in Eq. (|56|l . which are subsequently 
substituted in Eq. (|43|l to calculate Re £. The result is 



ReC 



c yrr 



1 + 36- 2 / 



Sirae- 2 / f 2 e 

n 2 -a 2 Q 



2 c -2 



(60) 

Thus the influence of the field is given by a sum of 
three contributions which can be extracted by consider- 
ing the frequency and field dependences of the response. 
All these contributions come about as small corrections 
~ e~ 2 , although the field itself (near the surface) is not 
small. This is because the Meissner layer is very thin. 



V. CONCLUSIONS 

I have considered the microwave response properties 
of a slab made of a material with coexisting super- 
conductivity and magnetism and calculated its surface 
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impedance for the parallel and perpendicular geometry. 
The impedance, and especially its real part, shows a rich 
behavior around threshold frequencies, related to spin- 
wave propagation. In the stiff regime, the propagation 
threshold is wj, above which Re£ shows a square- root 
singularity. On the other hand, in the soft regime this 
behavior is observed near a lower threshold u) m , while 
around u>f, Re£ has a square-root singularity below it 
and a linear behavior above it. In addition, there is a res- 
onance peak at Uh in both regimes, related to excitation 
of a surface spin wave^S The high-frequency response is 
mostly imaginary, being dominated by the superconduc- 
tivity, with a small real correction ~ oj~ 3 / 2 . 

This behavior is common for both geometries in the 
field-free situation, when no magnetic induction exists at 
the sample surface. However, when such an induction 
does exist, its influence is completely different in these 
geometries. In the perpendicular geometry, the magnetic 
induction penetrates the slab in form of vortices. This 
leads to a change in the effective anisotropy constant, 
resulting in a uniform shift of all thresholds and allowing 
one to detect the resonances by sweeping the magnetic 
field at constant frequency. However, it also leads to 
an increased dissipation, smearing the singularities and 
eventually pushing the response into the metallic regime. 

On the other hand, in the parallel geometry, the mag- 
netic induction is screened inside the Meissner layer. We 
found the influence of the induction in the regimes of very 



low and very high magnetic stiffness. In the former limit 
the response near u> m || is strongly direction-dependent: 
for B and Mo in the opposite directions, the response 
is not changed significantly, while if they are parallel to 
each other, Re£ is exponentially suppressed due to a po- 
tential barrier formed in the Meissner layer. In addition, 
it also shows a singularity ~ [lo — uj h y (0)] -2 , once the 
induction is strong enough so that ojh \\ (0) > ui m \\ ■ In 
the opposite very stiff regime dominated by strong su- 
perconductivity, the induction is screened out very fast 
and its influence is weak. Interestingly, the directional 
dependence of the response is opposite to the previous 
case, so that Re((B > 0) > Rc((B < 0). 

These results demonstrate that the magnetic and su- 
perconducting properties of SCFM's (including uncon- 
ventional superconductors with orbital magnetism) can 
be effectively investigated by microwave response mea- 
surements and provide a theoretical basis for such exper- 
iments. 
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